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Abstract 



We present a detailed study of the transport properties of a superconducting 
Pb/Cu microdot with a 2 x 2 antidot cluster. The superconducting-normal 
(S/N) phase boundary, critical currents and current- volt age characteristics of 
this structure have been measured. The S/N phase boundary as a function of 
field B [Tc{B)) reveals an oscillatory structure caused by the limited number 
of possible vortex configurations which can be realized in these small clus- 
ters of pinning centres (antidots). We have analyzed the stability of these 
configurations and discussed the possible dissipation mechanisms using the 
critical current {Jc{B)) and voltage-current (y{I)) characteristics data. A 
comparison of the experimental data of Tc{B) and Jc{B) with calculations 
in the London limit of the Ginzburg-Landau theory confirms that vortices 
can indeed be pinned by the antidots forming a cluster and that the ground- 
state configurations of the vortices are noticeably modified by sending current 
through the structure. The possibility of generating phase-slips as well as mo- 
tion of the vortices in the 2x2 antidot cluster has also been discussed. 

Pacs: 74.25.Dw, 73.23-b, 74.76-w 
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I. INTRODUCTION 



Laterally modulated superconductors have been intensively studied during the past few 
years (for reviews see |jl|J3). The interest towards these materials was stimulated by the 
advances in the microfabrication techniques which enabled the production of sub-micron 
structures with the relevant dimensions of the order of the superconducting coherence length, 
^(T) and penetration depth A(T). These microstructures impose certain constraints on the 
behaviour of the superconducting order parameter by confining the superconducting con- 
densate within the sample's boundaries, including those introduced through lateral nanos- 
tructuring. Superconducting wire networks and Josephson junction arrays ( JJA's) are some 
examples of the laterally modulated films, where the applied magnetic field provides a con- 
tinuously changing length scale which can explore commensurability between the lattice of 
vortices and the underlying artificially introduced lateral lattice [^]-|^. 

In networks, fluxoid quantization governs the properties and defines the allowed quantum 
states. Circular currents (vortices) are induced around the cells forming the structure, thus 
defining sets of specific vortex configurations. 

Recently, a new class of superconducting systems was introduced, consisting of a regular 
pinning array formed by a lattice of antidots (i.e. submicron holes) made in a type-II 
superconducting film At temperatures close to Tc superconductors with an antidot 

lattice behave like weakly coupled wire networks P], where vortex depinning dominates 
the dissipation Well below the transition temperature, where the fluxoid quantization 

condition implies flux quantization at the antidots, new features appear in the magnetization 
@J^, critical currents ||T^ and S/N phase boundaries 0] which have been attributed to the 
formation of multiquanta vortices at the antidots. 

Taking into account the vortex-vortex interactions in samples with ~ 10^ antidots is, 
however, not a simple exercise due to the very large number of the interacting vortices. 
From this point of view, a microdot with an antidot cluster (2 x 2, 3 x 3, etc.) with a 
small number of interacting vortices is a very promising "intermediate" system between a 
single superconducting loop with a finite strip width and a superconducting film with a huge 
array of antidots. At the same time, using finite microdots with four antidots (the antidot 
cluster), we still keep the most essential features of the vortex- vortex interactions in presence 
of antidots. The reduced number of interacting vortices simplifies exact calculations which 
can also be extrapolated for the analysis of the vortex behaviour in substantially larger 
antidot arrays. 

In this paper, we study the transport properties of such a model superconducting mi- 
crosquare containing four antidots. Further on, for simplicity, we shall call this structure 
a "2 X 2 antidot cluster". This system can be considered as four unit cells of the regular 
square lattice of antidots and the initial simplicity makes the 2x2 antidot cluster a good 
candidate for achieving an insight into the stable vortex configurations expected for larger 
systems. Only a very limited number of vortex configurations are allowed for the 2x2 an- 
tidot cluster and therefore both the experimental results are less difficult to be interpreted 
and also a computational analysis is easier to perform. Additionally, the limited number of 
vortex configurations makes these structures also interesting for flux logic applications Jill ■ 

In this paper we focus on measurements of the superconducting/normal (S/N) phase 
boundary, the magnetoresistance, the critical currents and the V{I) characteristics of a 
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superconducting microsquare with four antidots. These experimental results have been 
compared with calculations in the London limit of the Ginzburg-Landau theory and in 
the framework of the de Gennes-Alexander model EM. The comparison has demonstrated 



that several well defined vortex configurations can be induced in these samples by properly 
tuning the magnetic field and the temperature. At very low transport currents, only four 
of these configurations are stable. By increasing the transport current, however, other flux 
phases, which are unstable without transport current, can be realized. Based on the V{I) 
characteristics data, we also consider how the occupation of antidots by vortices influences 
the onset of dissipation. 

The paper is organized as follows. Section |I| describes the sample preparation, the ex- 
perimental techniques and characteristic superconducting parameters of the Pb/Cu 2x2 
antidot cluster. In section |111 A| , we present the results on the S/N phase boundary Tc{B) 
for the 2x2 antidot cluster and a reference superconducting microsquare without antidots. 
The calculation of the S/N phase boundary for the 2x2 antidot cluster in the London limit 
and in the framework of the de Gennes-Alexander model is described and compared 
with the experimental data. The discussion of possible effects of disorder is also presented. 
Section [1IIB| comprises the experimental results of the critical currents and V{I) character- 
istics. The V{I) data are used to determine the dissipation related to the vortex motion in 
the 2x2 antidot cluster. 



II. EXPERIMENTAL 

A. Sample preparation and characterization 

The 2x2 antidot clusters consisting of microsquares of 2x2 /im^ with four antidots (i.e. 
submicron holes of 0.53x0.53 /im^ and center to center distance of 1 /xm) were written by 
e-beam lithography in positive PMMA resist onto Si02 substrates. Afterwards, a bilayer 
consisting of 50 nm thick Pb and 17 nm thick Cu was evaporated in one single run in a MBE 
apparatus at pressures of 5xl0~^ Torr. The Si02 substrates were N2 cooled during evapo- 
ration which decreased the film roughness down to 1.4 nm and suppressed the interdiffusion 



at the Pb/Cu interface [|T2|. After the lift-off processing, the samples were characterized by 
X-ray, AFM and SEM. The X-ray results show that the Pb/Cu bilayers are polycrystalline 
with a preferential growth of the Pb and Cu in the (111) direction. An AFM picture of the 
2x2 antidot cluster of Pb/Cu with four leads attached for electrical connection, is shown in 
Fig. |l|a. Also shown (Fig. |I]b) is the AFM image of a reference sample which consists of a 
Pb(50 nm)/Cu(17 nm) microsquare of 2x2 /im^ without antidots. 

The top Cu layer of 17 nm, evaporated onto the Pb film, was used to protect the Pb 
film from oxidation and to enable electrical connection to the experimental apparatus using 
a wire bonding technique through the 150x150 /im^ electrical pads of the sample. Besides 
that, the Cu layer changes also the coherence length and the penetration depth of the Pb 
layer. We discuss the influence of this Cu layer on the superconducting properties of the 
Pb/Cu bilayers in terms of proximity induced superconductivity in the next section. 

The electrical transport properties were measured in a '^He cryostat using the four probe 
technique. The magnetic field was applied perpendicular to the film surface and the temper- 
ature stabilization was better than 0.4 mK. Since the mesoscopic samples are very sensitive 
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and easily destroyed by electrostatic charges, all necessary grounding precautions were taken 
and 1.2 kQ resistors were connected in series with each lead to damp parasitic voltage peaks. 
The resistance measurements where performed using an ac-resistance bridge whereas the 
V{I) characteristics and critical current measurements where carried out using a dc-power 
supply and a nanovoltmeter. The superconducting/normal (S/N) phase boundary was de- 
termined automatically by keeping the sample resistance at a specific fixed value (usually 
at 10 % of the normal state resistance) and varying the magnetic field and temperature. 
Similarly, the critical current measurements as function of field were performed by keeping 
the sample at a certain fixed voltage and varying the applied current and magnetic field. 



B. Superconducting parameters of the Pb/Cu bilayer 

In order to study the Pb(50 nm)/Cu(17 nm) bilayer which was used for the 2x2 antidot 
clusters, we prepared also a single Pb(50 nm) and Cu(17 nm) film and measured their 
properties. Below, we present the summary of these results. 

Comparing the resistivity due to electron-phonon scattering, for the Pb(50 nm) film, 
P300K — Ptk = 27.8 jjiVLcm and for the 17 nm Cu thin films psoox — Ptk = 12.47 fiQcm, with 
the published bulk values of 21 jiVLcm and 1.7 fiQcm respectively ^M, we conclude that the 



resistivity of the Pb film approaches the bulk value whereas the resistivity of the Cu film is 
enhanced with respect to the bulk value due to the finite size effect [jl8|. The mean free paths 



at 7K determined from the resistivity values are Iph = 33 nm and Icu = 7.2 nm respectively. 
Thus, both single films, Pb (50 nm) and Cu (17 nm) are in the dirty limit since Ipb < 
and Xcu >> 1 JT^I respectively, where = 83 nm is the BCS coherence length of the Pb 
||T6| and Xcu = i^vp) / {^nTlcukB) = 38 at 7 K is a dimensionless impurity parameter. 
The superconducting properties of the S layer are characterized by two important lengths, 

the temperature dependent coherence length, ^s'(T) = 0.85-^^^= and the penetration 
depth, Xs{T) =0.66Al ^^^^^^—^ in the dirty limit. Here is the London penetration depth 

y/l-T/Tc 



|16| . For the Pb (50 nm) film at T= K, these two quantities take the values ^pb{0) = 45 
nm and Xpb{0) = 39 nm, respectively. Experimentally, we can also determine ^pb{0) from 
the relation between the perpendicular second critical field Bc2,s in a S layer and its in-plane 
superconducting coherence length, ^5(0), given by 

5c2,s(0) = ^ f ° (2.1) 
2<5(0)^ 

where $0 = h/2e is the fiux quantum. The coherence length obtained in this way 
is .^p;,(0)=36 nm. Now, we can determine the Ginzburg-Landau parameter upb = 
Xpb{0) / C,pb{0) = 1.08 and since upb > 0.7 we conclude that the single Pb layer is a type-II 
superconductor. 

By covering the Pb with a Cu layer, a superconducting Pb/Cu bilayer is obtained with 
different characteristic parameters than those of the single Pb film. We therefore define the 
effective superconducting parameters as the ones that can be determined from measurements 
on a Pb/Cu bilayer. 
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The Pb/Cu antidot clusters had a superconducting transition temperature, Tcj^g = 
6.05K, whereas the reference sample (i.e. square of 2x2 nm^ without antidots) had a 
^cjvs = 5.55K. The systematic difference in critical temperature between the perforated and 
reference sample is probably related to a difference of the electrical properties of the interface 
between the Pb and Cu layers. The possibility for an oxidation at the Pb/Cu interface is 
higher in perforated samples. This results in a smaller proximity coupling and, possibly, a 
higher effective critical temperature. Nevertheless, it is obvious that in both cases the top 



Cu layer decreases T^j^^ below the = 7.2K of the bulk Pb due to the proximity effect |]T7 . 

From the measured critical fields of the Pb(50 nm)/Cu(17 nm) bilayer, we determined 
an effective superconducting coherence length as ^pf,/CM(0)=39 nm, using Eq. pTT] . 

One may expect that an effective penetration depth for the (NS) bilayer, Xns{0) could 
also be defined. Especially, since there should exist a penetration depth, AAr(O), determining 
the decay of the magnetic field in the normal layer with proximity coupling, which should 
be different from A5(0). The exact definition of these two quantities, AAr(O) and Xns{^), is 
not obvious. Since we know that \ns{^) > ^5(0), we will estimate a lower limit for AAr5(0) 
as Xpb{0) = 39 nm. Further in the paper (see section [lllB 2| ), a more accurate estimate for 
AAr5(0) will be given. 

With this information, we can determine a lower limit for the effective Ginzburg-Landau 
parameter Kpb/cu = Xpb/cu{0) /^Pb/Cu{0) as 1 for the Pb(50 nm)/Cu(17 nm) bilayer. Since 
i^Pb/Cu > 0.7 the Pb/Cu bilayer is a type-II superconductor. 

As shown above, by covering Pb with a thin Cu layer, we are not only protecting the 
former against oxidation, but we are also changing, in certain limits, the effective ^ and A. 

Pb/Cu bilayers are good candidates for studying proximity coupling effects [|1^] between a 
superconductor (S) and a normal metal (N) since interdiffusion hardly exists at the interface 
of these materials, especially when they are evaporated at low temperature. Therefore, each 
material remains within its own boundaries. 

The proximity effect []T^ refers to the induction of superconductivity in a contacting 
normal metal by the decay of the superconducting order parameter. A, from a supercon- 
ductor into a normal metal. Not only the order parameter in the superconductor decays 
when approaching the S/N interface from the superconducting size but also a non vanish- 
ing superconducting order parameter nucleates in the N layer close to the interface. Thus, 
one of the most important parameters defining the proximity effect between two layers is 
the so-called coherence length of the normal layer, ^n, which determines the length in the 
normal metal over which Cooper pairs can diffuse (see Fig 0). In the dirty limit and for our 
Cu parameters, ^cu = {hvplcu/ ^T^kBT)^^'^ = 28 nm at 6 K, where vp is the Fermi velocity 
of Cu. Note that ^cu > dcu where dcu = 17 nm is the Cu thickness, which means that 
the superconducting order parameter is finite over the full thickness of the Cu layer in the 
temperature interval at which experiments on the Pb/Cu antidot cluster were performed (5 
K < T < 6 K). 

Another very interesting parameter in proximity induced superconductivity is the ex- 
trapolation length [|13[, b = ■yC,cu(T) coth.{dcu/ ^(T)) = 98 nm at 6 K in the dirty limit (see 
Fig. where 7 = pcu/ PPb- This parameter 6 is a measure of the unfavourable influence of 
the N layer on the superconductivity of the S layer. If b << C,s{T) the effects of proximity 
coupling are important whereas if 6 >> ^s(T) the superconducting order parameter at the 
interface is almost not changed and the effects of proximity coupling are negligible. In our 
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case at 6K, 6 = 98 nm, $,pb{6K) = 108 nm, thus the proximity effect has an intermediate 
strength in the Pb(50 nm)/Cu(17 nm) bilayer. 



We can classify our samples as being in the Cooper limit |T3[ since they fulfill the con- 
straints dcu < icu and dpi, < ^o- In this limit, A(x) can be taken as a constant over the 
individual layers and discontinuous at the interface. The perpendicular second critical field, 
Bc2,NS for a (NS) bilayer of thickness [d^ + ds) can then be determined from the second 
critical field of the single superconducting layer, -Bc2,s, by the dirty limit expression WW^, 



B....s{0) = , f (22) 

1 + rjdN/ds 

where rj = ps/pN for a specular scattering at the interface. This expression is in good agree- 
ment with our experimental results on bilayers with different Cu thicknesses. In particular, 
for the bilayer Pb(50 nm)/Cu(17 nm) used for the fabrication of the antidot cluster, Eq. |2.2| 
reduces to: Bc2,ns{0) = tt-Bc2,5(0) with a=0.91. Experimentally, we have obtained q;=0.86 
which is in agreement with the calculated value within a 5%. Taking into account that the 
scattering at the interface is probably not completely specular, the obtained agreement is 
quite good. 

Combining Eq. |2.1| and Eq. p.2| , one can obtain the ratio between the superconducting 
coherence length for the Pb(50 nm) layer and that of the Pb(50 nm)/Cu(17 nm) bilayer, 
Cpb/Cu{0) / ^Pb{0)= 1-05 . Using the value ^5(0)= 36 nm for the Pb(50 nm) film, we expect, 
according to this relation, that ^pb/cu = 1-05 ^^(O) = 38 nm for the Pb(50 nm)/Cu(17 nm) 
bilayer film, which is within 5% agreement with the 'experimental' value determined from 
the measured second critical fields of bilayers of that specific thickness, ^pb/cu{0)=39 nm . 

In summary, from the analysis of the superconducting properties of the Pb/Cu bilayer 
we can conclude that: i) The Pb(50 nm)/Cu(17 nm) bilayer shows proximity induced super- 
conductivity of an intermediate strength, ii) The proximity induced superconductivity pen- 
etrates through the whole thickness of the Cu layer, iii) The effective parameters, -Bc2,Afs(0), 
^Ns{^) arid Ajv5(0) for the bilayer films have been determined from the experimental results. 



III. RESULTS AND DISCUSSION 
A. The S/N phase boundary 

1. Experimental results 

Fig. Ishows the (S/N) phase boundary, AT^{B)=T^{0) -T^{B), of the Pb(50 nm)/Cu(17 
nm) 2x2 antidot cluster measured with the criterium of 10% of the normal state resistance. 
Pronounced periodic oscillations of Tf.{B) are observed every 26 G. Defining a flux quantum 
per antidot as $0 = h/2e = B ■ S, where S is an effective area per antidot {S = 0.8/im^), 
the oscillations with periodicity of 26 G can be correlated with a magnetic flux per antidot, 
$ = n^o, where n is an integer number. In each of these 26 G periods, smaller dips appear 
at approximately 7.5 G, 13 G and 18 G. The smaller dips correspond to approximately 
$/$o = 0.3, 0.5 and 0.7. 

Superimposed with these oscillations, a parabolic background is observed, reflecting the 
B{T) dependence of the second critical field of the quasi-one-dimensional stripes constituting 
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the antidot cluster, = (V^^o) / (27rW(^(T)) |T9| where w is the width of the stripe. From 
this formula, an effective coherence length C,ns{^)= 35 nm is obtained for this particular 
sample with Tc =6.05 K. Above 52 G, deviations from the main 26 G periodicity and 
intermediate substructure are observed. 

In the inset of Fig. ^, two magnetoresistance curves, R{B), measured at Ti = 6.009 K 
and T2= 5.972 K are shown. At Ti only minima at nx 26 G are distinctly observed. At T2 
an additional substructure appears in each period. The magnetoresistance results confirm 
the features observed in the Tc{B). A cut of the Tc{B) boundary at high temperature (i.e. 
T ~ Ti) reveals only the main 26 G oscillations, whereas a cut at lower temperatures results 
in the full substructure. 

Fig. ^ shows the (S/N) phase boundary, ATc{B)=Tc{0) — Tc{B), of the reference sample 
(i.e. a superconducting microsquare without antidots (see Fig. |I|b)) measured at 10% of the 
normal state resistance. In this case, oscillations are also observed however they are clearly 
different from those presented in Fig. 0. The first oscillation of ABi ~ 16 G is followed 
by quasi periodic oscillations with decreasing AB {AB2 = 9.9 G, AB3 = 8.2 G, AB4 = 
7.5 G, AB = 7.1 G) In addition, the superimposed background of Fig. ^ follows a linear 
relation instead of the parabolic dependence seen in Fig. |^. The oscillations are caused by 
the confinement of the superconducting condensate in the dot as revealed by other studies 
pO| ||39|| . The remark must be made that the first two periods are noticably larger than 



predicted in Ref. |p9|.The other periods correspond to the calculations within 15% The 
linear background observed in Fig. ^ is related to the second critical field of a 2D system, 
B"^^, is given by the expression B"^^ = $o/2vr,^(T)^. From this linear background an effective 
coherence length, ^Ar5(0)= 41 nm is determined for this specific sample with Tc= 5.55 K. 

In the inset of Fig. |, the magnetoresistance curves for the reference sample are shown at 
four different temperatures. For any of these R{B) curves, two large drops of resistance are 
observed when decreasing the magnetic field. In addition, for each curve and at low fields, 
one single small dip appears. This dip shifts to higher magnetic fields when decreasing the 
temperature and is responsible for the oscillations observed in Tc{B). On the contrary, the 
two drops of resistance may be explained by the influence of the four narrow superconducting 
leads which are attached to the microsquare (two for sending the transport current and two 
to measure the voltage). These leads are quasi-one-dimensional wires and therefore they 
become superconducting below the field BI2 = (vT2$o)/27rw^(T). However, the dot itself 
has a second critical field transition given by -B^^ = <l'o/(2vr^(T)^), as mentioned above. 

Taking into account that the transition to the superconducting state takes place at higher 
magnetic fields for the ID system (i.e. leads) than for the 2D system (i.e. microsquare), 
one expects to attain the following situation: when decreasing the magnetic field, starting 
from the normal state, a first resistive drop is observed, caused by a transition to the 
superconducting state of the narrow current and voltage leads. The fact that the resistance of 
the microsquare is affected by this transition, does not mean that the "four point" resistance 
measurement contains a portion of the leads, but arises from those parts of the cluster which 
become superconducting due to their proximity to the leads. The theoretical critical field 
values for the leads, as calculated from the the above formula for Bl^ using the effective 
coherence length C,ns obtained from the linear background of Tc{B), are marked by the open 
squares in the inset of Fig. ^j. The first resistance drop of the dot takes place at slightly 
lower fields than the N/S transition of the leads, since they are only in proximity coupling 
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with the dot. Therefore, the squares indicating B^^ fall on a horizontal straight line above 
the onset of the resistance drop, and not in the middle of the transition. 

Further decreasing the magnetic field, we cross the line B"^^ at which the microsquare 
becomes superconducting and a second resistance drop is initiated (see open circles in the 
inset of Fig. ^). In this part of the curve, small dips appear in the magnetoresistance due 
to the dot geometry. These dips are responsible for the oscillations observed in the Tc{B) 
curve. Similar two-step R{B) transitions have been observed by Chi etal. who studied 



narrow ID superconducting wires connected to large 2D contact pads. In their case, the 
proximity coupling between the ID wires with a large critical field and the 2D pads with a 
small critical field gave rise to a two-step behaviour comparable to our results. 

We have to note that the expression used to calculate the B'^2 values corresponds only 
to the linear Tc{B) background; the periodic oscillations superimposed with it are not con- 
sidered in this expression and thus, it is not surprising that the calculated values for B^f* 
(see open circles in the inset of Fig. ^ do not lie on one horizontal line. However, if we 
determine the critical magnetic field values from the measurements of Tc{B) (see Fig. 4) 
where both the background and oscillations due to the dot geometry are considered at a 
criterium of 10%, we obtain the open triangles shown in the inset. Their position agrees 
with all the expectations, i.e. the triangles lie on a horizontal line which is at 10% of the 
normal state resistance. 



2. Calculations of the Tc(B) phase boundary 

In order to calculate the S/N phase boundary of the 2x2 antidot cluster, we have approx- 
imated the cluster geometry by a square network of quasi-one-dimensional superconducting 
wires with a width w and a length / (see Fig. H). Adjacent nodes are coupled by a 
supercurrent depending on the gauge-invariant phase difference between the nodes : 

7i,- = 0j + ^ TA-dl (3.1) 

where is the phase of the superconducting order parameter at site i, $o = h/2e the flux 
quantum, A the magnetic vector potential and d\ a segment along the wire. First of all, 
we have used the approximation developed for weakly coupled wire networks |^J^J2^|TI|] , the 
"interacting loop-current (ICL) model", which neglects the variation of the order parameter 
along the superconducting wires of the network and only considers the phase variations. This 
simplification is therefore equivalent to the London limit of the Ginzburg-Landau theory, 
and leads to a linear current-phase relation : 

<^owd 

where lij is the supercurrent from i to j, dis the wire thickness and A is the superconducting 
penetration depth of the material. The resulting kinetic energy for each wire can be 
expressed as : 
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and is a quadratic function of phase differences. The prefactor in Eq. ^]3] defines the couphng 
strength between the adjacent nodes and is denoted as : 

To obtain the possible energy states of the cluster, all phase differences •jij have to be 
determined. This is done by imposing current conservation at each node : 

Y.h, = ir (3.5) 

j 

where If ^* is the transport current fed into node (i) , and by applying the fluxoid quantization 
condition to each individual loop : 

$ 

^7ij = 27r(nfc - — ) (3.6) 

where the sum is taken over the wires of the fc*^ cell, (A; = 1. . . 4) is the flux quantum 
number of the k^^ cell, and $ the magnetic flux. An estimation of the self-inductance 
for an individual cell of the antidot cluster gives C ~2 pH. Since most of the results are 
obtained close to Tc, the flux generated by the circulating currents is negligibly small. As 
a consequence we have not taken the self- and mutual inductance effects into account in 
our theoretical analysis and $ is considered to be equal to the external magnetic flux. By 
solving Eqs. |3.5| and p.6| for the 7jj's at a given magnetic field and summing the corresponding 
energies Eij over all branches («, j) of the system, we obtain the total kinetic energy of the 
2x2 antidot cluster : 

E = \T. Ej^l (3.7) 



Since the phase differences are found from Eqs. p.5| and the total energy E for each set 
of quantum numbers Uk (fc=l,. . . ,4) becomes dependent on the magnetic flux through the 
cells. In Fig. |]a the total energy E is plotted as a function of the magnetic flux for the 
case of /f^* = and for all wires with identical lengths and widths. We have omitted the 
parabolic contribution arising from the finite width of the wires. Since the energy is periodic 
in $ with a period $o, we have only plotted the first period. Each set of {k = 1. . .4) 
values corresponds to a different vortex configuration and gives rise to an energy branch 
which is quadratic with the magnetic flux. Due to the symmetry of the structure, most of 
the branches are degenerate and only six principal parabolae can be distinguished. The S/N 
phase boundary is given by the branch which has the lowest energy for a given value of the 
flux 1T§: 



where V is the total volume of the structure, Tc is the critical temperature at zero field 
and ^(0) is the coherence length at T = 0. The phase boundary is therefore composed 
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of five branclies and has three minima at = 0, 0.5, 1 and four cusps at $/$o = 

0.3,0.37,0.63,0.7 (see Fig. |a). At the first branch, around $ = 0, all the = {k 
= 1...4) and a screening current is present only at the edge of the structure. When $ 
increases, the fourfold degenerate state with one vortex at one of the cells becomes more 
stable. The transition occurs at ^/^o = 0.3. At $/$o = 0.37 a doubly degenerate diagonal 
state becomes favorable, where two vortices occupy one of the diagonals of the cell. Note 
that this state (which is the analog of the "checkerboard" vortex configuration in antidot 
lattices p3[) has a lower energy than the configuration where two adjacent cells are occupied 
by vortices (parallel state). The diagonal state could be of interest for flux quantum logic 
applications In higher fields, a third vortex enters the cluster at $/$o = 0.7 and finally 
around $/$o = 1 all the cells are filled and only an edge current is flowing now counter 
clockwise. 

It is interesting to note that Eqs. |3.2| to |3.7| are similar to the equations describing the 
supercurrent and the energy in Josephson junction arrays [Q. In JJA's however, the current- 
phase relation is sinusoidal and only the phase differences over the Josephson junctions have 
to be considered. It turns out that the differences between the linear and the " J JA" models 
are not very important when the static properties in a magnetic field are considered. For 
comparison we show in Fig. ^ the E{^) curve for a JJA together with the one obtained for 
the linear relation. Nearly identical energy branches are indeed found caused by the obvious 
fact that the sinusoidal current-phase relation can be approximated by a linear one, for small 
phase differences. 

For completeness, we also show in Fig. |^ the curve predicted by the linearized de Gennes- 
Alexander formalism for strong-coupling one-dimensional wire networks, which incorpo- 
rates phase and amplitude variations of the order parameter along the wires of the network. 
This formalism assumes that the diameter of the superconducting strands which constitute 
the network is smaller than both the coherence length and penetration depth. In that way, 
the problem becomes 1 D, and the order parameter at a position s on a wire connecting 
two adjacent nodes a and b (see Fig. ^ can be expressed as follows: 



sm 



\I^aSin 



+ e '^"^ sin( 



(3.9) 



where \l/a and \l/f, are the order parameters at the nodes a and b, lab the distance between 
points a and b and '^ab the line integral of the vector potential multiplied by 27r/$o- From 
the second GL equation it is possible to determine the supercurrent density through the 
branch. It depends sinusoidally on the phase difference between the nodes a and b, similar 
to the current through a Josephson junction: 



J, 



ab 



hTr I UTr I S^T^{<Pa - (pb - lab) 
\^n\ \^b\ 



sm 



(3.10) 



Applying the standard boundary conditions at a node a, gives rise to the Alexander node 
equations: 



E 



cot 



/ 

''an 



sm 







(3.11) 
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which express the current conservation at node a similar to Kirchhoff's current law for 
resistor networks. The sum is taken over all nodes n which are nearest neighbours of node 
a (see Fig. |^). If this equation is written for all nodes in the network, this leads to a 
characteristic determinant which must be zero for the existence of a non-trivial solution. 
We obtain 



1 + 2 cos — 3 cos —j-=;r ■ 1 + 2 COS — h 3 COS —j-=;r ■ (3.12) 



2U . 27r<l>\ / 2U . 




where $ is the flux of the applied magnetic field per antidot. Equating the first factor to zero, 
one obtaines a curve defining the first and the fifth branch of the Tc{B) phase boundary. 
Likewise, the middle three branches are determined by the other three factors (see Fig.|^). 

This linearized approach is valid only close to Tc. At lower temperatures, the non-linear 
version of this model should be applied pSj^ . 



The qualitative agreement between the ILC model (see Fig. ^(a)) and the experimental 
curves of Figs. 1(6) and (c), where the parabolic background has been subtracted, is fairly 
good. The measurement of curve (6) was performed with an ac current of 1 /lA, the one 
in (c) with 3 fiA. The five parabolae corresponding to the different states are found back 
in the experimental plots. An important difference is however, that two maxima appear 
in the experimental phase boundaries at $/$o — 0.2 and 0.8 which are not reproduced in 
the calculations. The ATc($) dependence for the states with all = or all = 1 thus 
seems steeper than predicted by the theoretical model and it also has a larger amplitude. 
Strangely, the agreement with the theoretical curve, which is calculated for zero transport 
current, is better in curve (c), where the highest transport current was used. 

In addition, the experimental ATc($) (see Fig. ^ shows a considerable parabolic back- 
ground, some distortions and a disappearance of the substructure above 2$o- 

As it was mentioned previously, the parabolic background ("diamagnetic shift") is due 
to the penetration of the magnetic field in the volume of the wires and can be taken 

into account by adding the term : 

.E^.^^ (3.13) 

to the right hand side of Eq. |3.7| . The effect of this background is shown in Fig. ||. Curve (a) 
shows ATc(<l') obtained from Eq. |3.7| and Eq. |3.8| for the antidot cluster. Curve (6) considers 
the case of a finite linewidth (Eq. p.l3| ) with w = 0.35 fim for wires at the edges of the 
network and w = OA fim for the wires connecting the edges with the central node (see 
Fig. 1^). Note in curve (6) that as the field is increased, the ATc amplitude due to fiuxoid 
quantization in the structure becomes quite small with respect to the shift caused by the 
background and the amplitude of the oscillations gradually decreases. 

After discussing in detail the first period of the phase boundary Tc($/$o), "we would like 
to analyze briefly the possible reasons for the shift of the Tc($/$o) minima in the higher 
periods of the S/N phase boundary, i.e. for 1< $/$o <2, etc. (Fig. H). Taking into account 



11 



the "softness" of the current loops in a real structure, which, strictly speaking, is not a one 
dimensional network, the presence of disorder of these loops is quite probable. Therefore, 
one may assume that one of the reasons of the variation of the Tc($/$o) peaks from period 
to period can be related to an areal disorder of "soft" current loops in the 2x2 antidot 
cluster. 

Another important factor is the disorder arising from width inhomogeneity, structural 
defects, non identical electrical properties of each wire, which can lead to a modification 
of the critical current of the individual wires and to a distribution of the effective area per 
cell. The importance of areal disorder was already stressed in studies of the phase boundary 
in wire networks [0,0 where it was shown that areal disorder can lead to a decay of the 
oscillation amplitudes and in some cases even to beatings in the envelope. Measurements 
on Josephson junction clusters (JJC's) |^ showed that coupling disorder coming from an 



unavoidable spread in junction parameters is not averaged out as in large arrays, leaving a 
clear trace in the transport properties. 

Curve (c) in Fig. |] shows the influence of areal disorder on the theoretical phase bound- 
ary of the antidot cluster. The areal disorder was introduced by allowing the coordinates of 
the nodes to vary randomly within a circle of radius 0.1 fim around the node position for the 
ordered network. In this way a random distribution of lengths and cell areas is generated 
which lifts the degeneracy of the possible states and changes the relative positions of the 
different parabolic energy branches. For small fields (i.e. the first period), the deviations 
with respect to curve (6) are not very pronounced and it is still possible to identify all five 
parabolae forming the Tc{B) phase boundary. As the field increases, the oscillation ampli- 
tude gets smaller and the positions of the different branches shift, making the identification 
of the states less straightforward. From this analysis it is clear that disorder can indeed 
cause a shift of the peaks with increasing magnetic field, as it is observed experimentally. 
Finally, in curve {d) we have included the contributions of the parabolic background (curve 
(b)) and the areal disorder (curve (c)). A comparison of curve {d) with the experimental 
data of Fig. ^ shows that certainly these two contributions are playing a role in our measure- 
ments. Nevertheless, even by including these two effects in the model, it is still not possible 
to simulate the experimental curves completely. This is probably due to the fact that we 
have used a one-dimensional approach, which cannot take into account the 2D character of 
the structure. 



B. Critical currents and V{I) characteristics 

So far we have considered the effects of the vortex confinement by the 2x2 antidot cluster 
on the S/N phase boundary Tc{B). In order to demonstrate that the unique properties of 
the Pb/Cu 2x2 antidot cluster are not restricted to the S/N phase transition, we present 
below the critical current results and the V{I) characteristics measured at temperatures 
400 mK < Tc — T < 100 mK. We will show that at these temperatures the quantized states 
are still present and that the most stable ones at the S/N phase boundary do not always 
correspond to the states carrying the highest currents. In addition, we will demonstrate that 
a transport current j 7^ is able to lift some of the degeneracies of the vortex configurations 
at j = 0. The V{I) characteristics will be used to give a qualitative picture of the fiux line 
transport in these nanostructures in terms of phase-slip processes. 
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It is known that the dynamics of ID wires p3|,^ and arrays of ID wires are mostly 
governed by phase-shp processes when currents close to the depairing critical current are sent 
through the wires. Thus, one may also expect that signatures of these processes should also 
be present in the V{I) characteristics of the 2x2 antidot cluster. In such a phase-slip process 
the energy of the system is reduced by bringing a small spot of the ID wire (of the order 
of the quasiparticle diffusion length, Xq ^ ^{T)) momentarily to the normal state. During 
this process, the phase of the superconducting order parameter in that spot is changed by 
27r. To preserve the superconductivity in the sample in the presence of a large current, the 
phase-slip process repeats in time and the average period r between phase-slips is related 
to the voltage measured between the two ends of the wire through the Josephson relation, 
V = The points where the order parameter becomes zero and its phase shows jumps 

of 27r, are known as phase-slip centres (PSC). Obviously, if the wire already has some weak 
superconducting points, the PSC's will be localized at these spots. 

Several mechanisms of PSC formation have been reported. First, PSC's can nucleate 
when the current exceeds the depairing critical current of the ID wire, as described by the 
theory of Skocpol, Beasley and Tinkham (SBT) |3^. In this case the formation of 1,2,. . . n 
PSC's gives rise to a step-like V{I) characteristic. 

PSC's can also be formed as a result of thermodynamic fluctuations which take place 
with a probability proportional to exp {-5T/kBT), where 5T = wV2Ejl/3^{T) is the free 
energy barrier between the state before and after the phase-slip. The theoretical description 
of PSC's was developed by Langer and Ambegoakar (LA) [^. Taking into account the 



exponential decrease of the phase-slip probability with temperature, their model is only 
applicable in a very narrow temperature interval near Tc. 

A third mechanism has recently been reported by Giroud et al. for arrays of ID 
wires with localized vortices. They consider the possibility that when a current close to the 
depairing critical current is reached in one of the wires, the vortex feels a " Lorentz-like" 
force perpendicular to the transport current, which tends to move the vortex to the next 
cell where the process is repeated. This model assumes the nucleation of a phase-slip in a 
ID wire each time that a vortex crosses the ID wire. 



1. Critical currents: experimental results and comparison with the model 

The critical current density versus magnetic field, Jc{B), of the 2x2 antidot cluster was 
determined with a criterium of 3fiV as explained in section [11 A| . Then, at specific values of 



the magnetic field and temperature, the V{I) characteristics were measured to confirm the 
Jc(-B) results and study the vortex dynamics in the antidot cluster. 

Fig. ^ shows the Jc{B) curves measured at four temperatures. Note that clear maxima 
are observed at n$o and (?2 + 0.5)$o, and smaller inflections are detected around (n + 0.3)$o 
and (n + 0.7)$o (see Fig. |b and |c). Initially, the magnitude of the oscillations increases by 
decreasing the temperature (Fig. Pa,b) but they disappear almost completely when further 
decreasing the temperature (see Fig. ^i). A parabolic background caused by the magnetic 
field penetration in the elementary wires is again observed at all temperatures. 

The magnetic field values at which specific features appear in Jc{B) (Fig. ^)and Tc{B) 
(Fig. 0) are the same. Thus one may assume that the maxima observed in the Jc{B) of 
the 2x2 antidot cluster are related to certain stable vortex states revealed in the Tc{B) 
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oscillations. However the relative amplitude of these maxima in one period of $o is different 
for Jc(-B) than for Tc{B). In Tc{B) large minima were observed at ?2$o and three smaller 
minima of the same energy were obtained at around (n + 0.3)$05 (^ + 0.5)$o and (n + 0.7)$o- 
On the contrary, in Jc{B), the maxima at n$o are followed in magnitude by the maxima 
at (n + 0.5)$o whereas the features appearing at {n + 0.3)$o and {n + 0.7)$o have almost 
completely lost their amplitude. 

To estimate the field dependence of the critical current theoretically and compare it 
with the experimental observations we have used the model described in section P^II A2| and 
solved equations ^.5| and |3.6| for the case where the external transport current enters the 



structure at node (a) and leaves at node (6) (see Fig. The current-phase relation, given 
by Eq. |3.2| , has slightly been modified to take into account the order parameter depression 
in the strands when the current is close to the depairing current 1^^. We have used the 
relation : 

/.. = /o(^)(l-^(^f) (3.14) 

which is linear for small phase differences and becomes parabolic near = 7c, where 
7c = //(v^^) is the critical phase difference at which a phase-slip process occurs and 
/o7c = Ifep is the Ginzburg-Landau depairing current. Note that Eq. p.l4| is identical to 



the current-phase relation obtained for a long ideal weak link in the depairing limit ||30|| . 



The set of nonlinear equations for the unknown phase differences 7^^ obtained from 



Eqs. |3.5| and 3.6, was solved numerically using the standard Newton-Raphson method [31 



The critical current at a magnetic flux $ and a fixed state (fc = 1 . . . 4) was obtained 
as follows : first, the 'jij^s were initialized to zero, /^^* was fixed at a certain value and the 
corresponding phases were determined. Next, the external current was slowly ramped up 
and the same procedure was repeated until the current was too high to find any solution 
for the set of equations. The external current, above which no static solution exists, has 
been taken as the intrinsic critical current of the structure. The self-field induced by the 
transport current was estimated to be not higher than 0.1 % and was therefore neglected in 
the analysis. 

Curve (a) in Fig. |10| shows the theoretical Ic{^) obtained with Eq. p.l4| and 7c = vr/2. 
Only the states with the highest critical current at that particular <l> are shown and we 
have restricted ourselves to the first period, without considering any disorder. Because of 
the current injection at node (a) the symmetry is broken and the resulting vortex config- 
urations differ from the ground-state configurations at zero applied current, discussed in 
section [111 A 2| . 

The possible states are displayed schematically in Figure 0^a). Only two states with 
one vortex in the structure are possible, instead of the fourfold degenerate ground-state 
obtained for the case of I^^^ = 0. Near $ = |<l>o the supercurrent is carried by a state where 
the vortices occupy the second row of the cluster (parallel state) instead of being located on 
the diagonals (checkerboard configuration). 

The reason for the substantial modification of the stable vortex configurations is of course 
that the external current is added to the circular currents which flow to satisfy the fiuxoid 
quantization. The total amount of current that can be injected before the structure depairs 
is higher if the transport and circular shielding currents are subtracted. For the case of large 
weakly-coupled wire networks |^ and inductive JJA's it was already predicted that large 
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transport currents make the checkerboard ground-state unstable. This leads to a state of 
parallel vortex rows which can move coherently (driven vortex lattices) over the underlying 
periodic array. 

Curve (6) in Fig. |l^ shows the static calculation for a JJC having a sinusoidal current- 
phase relation for comparison. Besides some minor differences concerning the exact positions 
of minima and maxima (compare curves (a) and (6) in Fig [1^), the same states can be 
identified in the Ic{^) plot. Note however, that in this case the normalizing current Ico is 
the Josephson critical current and not the depairing current. 

If in the case of the antidot cluster (curve (a)), 7c > 7r/2, the overall shape of the calcu- 
lated Ic{^) curve slightly changes and the difference with the JJC curve is more important, 
although still the same states carry the largest supercurrent. 

If we compare the theoretical curve (a) with the first period of the experimental data 
shown without the background for T=5.808 K in Fig. |lO|(c), we note that the curves are 
qualitatively similar. The same five states are clearly present. The state at $/$o = 0.5 has 
a higher than the states around ^/^o = 0.3 and 0.7 though it does not reach the same 
current value as predicted in the calculation, except for temperatures below approximately 
5.715 K. 

Below this temperature (T=5.715 K), the peaks at n$o flatten and the amplitude of 
the modulations of the critical current decreases with decreasing temperature. It seems 
that the intrinsic Ic{^) behaviour of the structure is cut off by a superimposed parabolic 
/c($) background (see the full line in Fig. 1). The disappearance of the substructure with 
decreasing temperature was also observed in wire networks and attributed to an increase 
of the energy barrier required to cross the superconducting wires. 

In our case however, it looks more like a cut-off rather than a continuous decrease. 
Therefore, we believe that the effect could be caused by the propagation of heat generated 
in the current leads. Since the antidot structure can be considered as a kind of parallel 
circuit having a critical current which is a factor approximately 1.7 times higher than the 
depairing current ( of a single wire (see Fig. [iy)(a), the current lead is probably in a 
resistive state while the structure is still not. We measured the voltage between the nodes 
(a) and (6) (see Fig. |^) thus the voltages appearing over the current leads should not have 
influenced the results. However, at lower temperatures (when 1^ is rather high), the actual 
heating due to this current cannot be excluded. If such heat propagates towards the antidot 
cluster, it could trigger a transition to a dissipative state earlier than expected. In that case, 
the /c(^) dependence would show a "cut-off" governed by the parabolic critical current 
dependence of the current-leads, I^i^)- 

Besides this influence of the current leads at high current levels, care should also be taken 
not to overestimate the limits of validity of the static simulation itself. Unlike the case of 
JJA's where the dynamics can relatively easy be taken into account by using the resistively 
shunted Josephson junction (RSJ) model in combination with the Josephson voltage relation 



3^, it is not so straightforward to describe the dynamics in wire networks where resistive 
PSC's nucleate and can be activated throughout the structure. The present model does 
not take into account any dynamic effects such as vortex motion or propagating PSC's and 
therefore, the critical currents obtained by the static simulation should be considered as an 
estimate of the upper limit, rather than an accurate determination. 
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2. V{I) characteristics: results and discussion 



Fig. O shows three V{I) curves measured at T=5.715 K for three magnetic fields as 



representative data. Note that in all the curves, several steps are observed which depend 
on the applied magnetic field. However, an enlargement of the lower part of the V{I) 
characteristics (see Fig. O) shows that a quasi-linear dissipative foot is clearly observed for 



the cases of $ 7^ ^'^'o- Note that this foot strongly depends on the magnetic field. In Fig. ^ 
we show three curves corresponding to $ = $0; ^ = 0.37$o and $ = 0.5$o for T = 5.808 K 
and T=5.715 K. Indicated by a horizontal dotted line is the 3/iV criterium used to determine 
the Jc($) of Fig. a 



When $ = n$0; no foot is observed at any temperature. Fig. |I3| shows the linear de- 
pendence of the critical current, I^^^, of the Pb/Cu 2x2 antidot cluster versus T/Tc for the 
measurements performed at zero magnetic field. Ic follows the (1 — T/TcY^"^ dependence ex- 
pected from the Ginzburg-Landau depairing critical current of ID wires. The extrapolation 
of the (1 - law to T/T^ = results in /c(0)=7 mA. 

On the other hand, an estimation of the theoretical Ginzburg-Landau value for the 
depairing critical current of a ID wire gives 



where w and d are the width and thickness of the wires respectively and ^(0) and A(0) were 
determined in section [II B[ As shown in the previous section (see Fig.^(a)), the depairing 



current of the structure is ll^{0) ~ 1.6/^^(0) = 72 mA, which is one order of magnitude 
higher than the experimental value obtained. It should however be noted that we have 
underestimated the value of A(0) by equating it to that of the single Pb layer, Xpb{0), while, 
in fact, we should use Xns > Xpb (see section [II B| ). From the above analysis, we can estimate 
an effective penetration depth for the Pb/Cu bilayer Xns{0) ~ 3Ap6(0). 

When $ 7^ $0) we believe that the field dependent linear foot which appears at low 
currents in the V{I) characteristics, is induced by phase-slip processes. More precisely 
by the mechanism reported by Giroud et al. |0] for arrays of ID wires which claims that 
vortices may move from one cell to the other in a direction perpendicular to the transport 
current in a kind of stationary "fiux-flow" regime. In this model a vortex jump from one 
antidot to another should lead to a phase-slip induced in the ID wire crossed by the vortex. 
Dissipation in arrays due to the vortex motion has also been reported by other authors . 
This mechanism may explain the disappearance of the resistive foot at $ = n$o- Since we 
have shown that at $ = ?t,<I>o an edge surface current is fiowing in the antidot cluster, no 
vortices are present inside the antidot cluster. Thus, if the dissipative foot is induced by the 
motion of vortices, no foot should be expected at $ = r2$o- 

The possibility of vortex motion was not considered for the critical current calculation in 
the previous section. There we determined the critical current of the structure using a static 
approach, i.e. the Ic was defined as the maximum I which could keep zero voltage through 
the structure, and still fulfill the fiuxoid quantization condition and current conservation 
equations. The real critical current values may therefore be lower than the ones obtained 
by means of the static approach. 



16 



The experimental evidence supporting the idea that motion of vortices in the antidot 
cluster induces phase-slip processes and thus gives rise to dissipation at / < Ic, is presented 
below. 



Fig. shows the differential resistance, dV/dl, as a function of the reduced flux, $/$o, 
determined at the onset of the dissipative foot in the V{I) characteristic for T=5.808 K and 
T=5.715 K (Fig. p^) . Since the V{I) characteristic for the dissipative foot is quasi-linear, 



we have extracted one single value for the differential resistance at each $/$o- Fig. WA 



shows the dV/dl for the $/$o interval where the current is flowing in the internal strands 
of the 2x2 antidot cluster. The interval where only an edge surface current is flowing 
in the antidot cluster has been omitted since there, no foot is observed. We have denoted 
by vertical dashed lines the $/$o values at which a change of vortex configuration was 
experimentally observed. Note that at both temperatures, the dV/dl can be defined by two 
levels: a low level which accounts for the one- vortex configuration (N=l) and three- vortex 
configuration (N=3) or also called one-antivortex configuration, and a high level given by 
the two- vortex configuration (N=2). This two level representation strongly suggests that 
the dissipative foot might be caused by a "flux- flow" motion of vortices as reported also by 
Giroud etal. for wire networks. In this model, a vortex crossing one cell is associated to 
a 2tt phase jump, which generates a voltage in a time r given by the Josephson relation, 

h2'n , , 

V = -- (3.16) 

If we assume that the stationary dissipation process can be described by a constant average 



vortex velocity v (as it the case of classical flux-flow in type-II superconductors ||3^ and 
JJA's [0), the voltage generated by the net motion of N vortices can be written as : 

V = A2^/ = W„!; (3.17) 

2e r a 

On the other hand, since this model associates a PSC to each vortex jump, we have equated 
the total average power V^(J — /c($)) generated in the cluster to the sum of powers dissipated 
by each individual PSC, Ppscs = NiQRn2XQ/ a, where here N is the number of PSC's, iq is 
the quasi-particle current flowing in each wire, i?„ is the normal state resistance of a single 
wire and Ag is the quasiparticle diffusion length. From this equality and assuming that 
~ (/ — /c)/3 |]3|,|33|, we have determined an average vortex velocity, 

Rn2Xri{I — Ir) 

^ 3.18 



which is proportional to the current (/ — Ic) ■ By substituting this velocity in Eq. p.l7| and 
calculating the differential resistance at the critical current (Jc), we obtain: 

2 NRnX*n 

= Q 3.19 

9 a 




From Eq. p.l9| and taking Ag ~ ^(T) and N=l for the one- vortex configuration and N=2 for 



the two-vortex configuration, we have determined the following differential resistances: At 
T=5.715 K, dV/dl ~0.1f2 for N=l and dV/dl ~0.2 f2 for N=2, which is in good agreement 
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with the values given in Fig. 0. At T=5.808 K, dV/dl ~0.12 Vt for N=l and dV/dl ~0.24 
for N=2. In this case, the values shown in Fig. |l4|are higher although their ratio still remains 
approximately correct. These calculations also point in the direction that the dissipative 
foot can be interpreted in terms of vortex motion. 

At higher currents, the V{I) characteristics do not show a continuous increase of voltage 
with increasing current but instead they are composed of several voltage steps (see Fig. 
The transition to the normal state of a ID wire is produced by the generation of steady PSC's 
which increase in number as the transport current grows. This dissipation phenomenon 
correlates the differential resistance with the number of PSC's created in the wire and 
therefore results in a V{I) curve with characteristic steps. Although the steps observed 
in the V{I) curve in Fig. |ll| are not steep like the ones reported for single ID wires ||33 



and infinite wire networks interpreted by the SBT theory, we believe that they are 
reminiscences of the resistive transition of the quasi-one dimensional wires constituting the 
antidot cluster. A non-strictly ID character of the wires forming the antidot cluster may be 
one of the reasons for the smoothness of the steps. 



IV. CONCLUSIONS 

We have studied the transport properties of a Pb/Cu 2x2 antidot cluster (a microsquare 
with four antidots) by measuring the superconducting/ normal phase boundary, critical cur- 
rents and V{I) characteristics. 

The Pb/Cu bilayer can be considered as a single superconducting entity with supercon- 
ducting parameters which are somewhat different than those of the single Pb layer. The 
presence of the antidots leads to a characteristic structure in the magnetoresistance and 
phase boundary caused by the formation of well defined vortex configurations. By compar- 
ison of the experimental results with calculations which approximate the cluster as a ID 
micro-network, we were able to identify the corresponding ground states. 

The formation of particular vortex states, as the magnetic field is varied, was also ob- 
served in the critical current versus field curves. Comparison with a static model showed 
that the current injection lifts the degeneracy and that the vortex states at larger current 
differ from the ground states, observed in the phase boundary Tc{B). 

The dissipative processes were probed for the different vortex occupations by means of 
V{I) measurements. We find clear evidence for a stationairy vortex motion at the onset of 
dissipation and the creation of phase slip centers at larger voltages and currents. 
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FIGURES 



FIG. 1. Atomic force microscopy picture of (o) the 2x2 antidot cluster and (6) the reference 
sample without antidots. 

FIG. 2. Schematic representation of the spatial dependence of the superconducting order pa- 
rameter at the interface between the superconductor (S) and the normal material (N) of the Pb/Cu 
bilayers studied in this paper. The two characterisic parameters for proximity induced supercon- 
ductivity, i.e. the coherence length in the normal layer, and the extrapolation lenght, b, are 
indicated. 

FIG. 3. Experimental S/N phase boundary of the 2x2 antidot cluster determined with a 
criterium of 10% of the normal state resistance. Two typical magnetoresistance curves are shown 
in the inset. The ac measuring current was l/zA. 

FIG. 4. Experimental S/N phase boundary of the reference sample determined with a criterium 
of 10% of the normal state resistance. The inset shows four magnetoresistance curves measured 
at T=5.503 K, 5.482 K, 5.462 K and 5.443 K. The ac measuring current was 1 /xA. The second 
critical magnetic fields for the ID transition (-Bc2^), the 2D transition, (B"^^), and the transition 
of the microsquare {B^^) are shown. 

FIG. 5. Schematic representation of the 2x2 antidot cluster with average dimensions indicated. 
The network approximation used in the one-dimensional model is drawn with a dashed line. The 
dots denote the nine nodes considered in the model. 

FIG. 6. (a) Theoretical phase boundary, AT($/<I>o), obtained in the London limit of the 
Ginzburg-Landau theory (solid line) . The total energy for the six vortex configurations is shown by 
a dot line. The dashed line indicates the non-stable 'parallel vortex configuration'. The schematic 
representation of the vortex states in each parabolic branch is also sketched by means of the cor- 
responding Uk quantum numbers. (6) First period of the experimental phase boundary shown in 
Fig. 1^ after subtraction of the parabolic background {lac = 1/iA). (c) idem as (6) but for lac = 3/iA. 

FIG. 7. Comparison of the theoritical phase boundaries obtained for the London (dashed line), 
Josephson junction array (solid line) and de Gennes- Alexander (dot line) approaches. 

FIG. 8. (a)Theoretical phase boundary for an ideal one-dimensional antidot cluster, [h) The- 
oretical phase boundary for an antidot cluster considering the effects coming from the magnetic 
field penetration in the quasi-one-dimensional wires, (c) Theoretical phase boundary for an antidot 
cluster taking into account the effects induced by areal disorder, (d) Theoretical phase boundary 
where the contributions shown in (6) and (c) are both considered. 
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FIG. 9. Critical current density of the 2x2 antidot cluster as function of the reduced flux, 
$/«>o, at four temperatures, (a) T=5.899 K, (6) T=5.808 K, (c) T=5.715 K and (d) T=5.656 K. 



FIG. 10. Critical current versus the reduced flux per cell for (a) the current-phase relation 
given by Eq. 3.14 and (b) for a sinusoidal current-phase relation. The vortex states generated in 
the 2x2 antidot cluster at high transport currents are shown, (c) Experimental critical current 
versus reduced flux per cell at T=5.808 K with the background subtracted. 



FIG. 11. V{I) characteristics of the 2x2 antidot cluster measured at T=5.715 K for zero 
applied magnetic field and a magnetic field of 3 G and 12 G. In all the curves similar characteristic 
steps can be distinguished. 



FIG. 12. Enlargement of the lower part of the V{I) characteristics for a reduced flux $ = $oj 
^> = 0.37$o and $ = 0.5$o at T=5.808 K and T=5.715 K. 



FIG. 13. Temperature dependence of the critical current measured at -B = showing the same 
behaviour as the depairing critical current of ID wires. 



FIG. 14. Differential resistance determined at the onset of the dissipative foot of the V{I) 
characteristics as a function of the reduced flux for T=5.808 K and T=5.715 K. The solid line is 
a guide for the eyes showing the two levels expected from the model of vortex motion induced by 
phase slips. The dot lines indicate the $/<&o values at which a change of vortex configuration is 
observed experimentally. 
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